Transcriptome of the dead: characterisation of immune genes and marker development from necropsy samples in a free-ranging marine mammal

Background Transcriptomes are powerful resources, providing a window on the expressed portion of the genome that can be generated rapidly and at low cost for virtually any organism. However, because many genes have tissue-specific expression patterns, developing a complete transcriptome usually requires a 'discovery pool' of individuals to be sacrificed in order to harvest mRNA from as many different types of tissue as possible. This hinders transcriptome development in large, charismatic and endangered species, many of which stand the most to gain from such approaches. To circumvent this problem in a model pinniped species, we 454 sequenced cDNA from testis, heart, spleen, intestine, kidney and lung tissues obtained from nine adult male Antarctic fur seals (Arctocephalus gazella) that died of natural causes at Bird Island, South Georgia. Results After applying stringent quality control criteria based on length and annotation, we obtained 12,397 contigs which, in combination with 454 data previously obtained from skin, gave a total of 23,096 unique contigs. Homology was found to 77.0% of dog (Canis lupus familiaris) transcripts, suggesting that the combined assembly represents a substantial proportion of this species' transcriptome. Moreover, only 0.5% of transcripts revealed sequence similarity to bacteria, implying minimal contamination, and the percentage of transcripts involved in cell death was low at 2.6%. Transcripts with immune-related annotations were almost five-fold enriched relative to skin and represented 13.2% of all spleen-specific contigs. By reference to the dog, we also identified transcripts revealing homology to five class I, ten class II and three class III genes of the Major Histocompatibility Complex and derived the putative genomic distribution of 17,121 contigs, 2,119 in silico mined microsatellites and 9,382 single nucleotide polymorphisms. Conclusions Our findings suggest that transcriptome development based on samples collected post mortem may greatly facilitate genomic studies, not only of marine mammals but also more generally of species that are of conservation concern

Tipping cycles

Instability in Jacobians is determined by the presence of an eigenvalue lying in the right half plane. The coefficients of the characteristic polynomial contain information related to the specific matrix elements that play a greater destabilising role. Yet the destabilising circuits, or cycles, constructed by multiplying these elements together, form only a subset of all the cycles comprising a given system. This paper looks at the destabilising cycles in three sign-restricted forms in terms of sets of the matrix elements to explore how sign structure affects how the elements contribute to instability. This leads to quite rich combinatorial structure among the destabilising cycle sets as set size grows within the coefficients of the characteristic polynomial

Transcriptomic SNP discovery for custom genotyping arrays: impacts of sequence data, SNP calling method and genotyping technology on the probability of validation success

Background Single nucleotide polymorphism (SNP) discovery is an important goal of many studies. However, the number of ‘putative’ SNPs discovered from a sequence resource may not provide a reliable indication of the number that will successfully validate with a given genotyping technology. For this it may be necessary to account for factors such as the method used for SNP discovery and the type of sequence data from which it originates, suitability of the SNP flanking sequences for probe design, and genomic context. To explore the relative importance of these and other factors, we used Illumina sequencing to augment an existing Roche 454 transcriptome assembly for the Antarctic fur seal (Arctocephalus gazella). We then mapped the raw Illumina reads to the new hybrid transcriptome using BWA and BOWTIE2 before calling SNPs with GATK. The resulting markers were pooled with two existing sets of SNPs called from the original 454 assembly using NEWBLER and SWAP454. Finally, we explored the extent to which SNPs discovered using these four methods overlapped and predicted the corresponding validation outcomes for both Illumina Infinium iSelect HD and Affymetrix Axiom arrays. Results Collating markers across all discovery methods resulted in a global list of 34,718 SNPs. However, concordance between the methods was surprisingly poor, with only 51.0 % of SNPs being discovered by more than one method and 13.5 % being called from both the 454 and Illumina datasets. Using a predictive modeling approach, we could also show that SNPs called from the Illumina data were on average more likely to successfully validate, as were SNPs called by more than one method. Above and beyond this pattern, predicted validation outcomes were also consistently better for Affymetrix Axiom arrays. Conclusions Our results suggest that focusing on SNPs called by more than one method could potentially improve validation outcomes. They also highlight possible differences between alternative genotyping technologies that could be explored in future studies of non-model organisms

Establishing RNAi in a non-model organism: The Antarctic nematode Panagrolaimus sp. DAW1

The Antarctic nematode Panagrolaimus sp. DAW1 is one of the only organisms known to survive extensive intracellular freezing throughout its tissues. Although the physiological mechanisms of this extreme adaptation are partly understood, the molecular mechanisms remain largely unknown. RNAi is a method that allows the examination of gene function in a direct, targeted manner, by knocking out specific mRNAs and revealing the effects on the phenotype. In this study we have explored the viability of RNAi in Panagrolaimus sp. DAW1. In the first trial, nematodes were fed E. coli expressing Panagrolaimus sp. DAW1 dsRNA of the embryonic lethal genes rps-2 and dhc, and the blister gene duox. Pd-rps-2(RNAi)-treated nematodes showed a significant decrease in larval hatching. However, qPCR showed no significant decrease in the amount of rps-2 mRNA in Pd-rps-2(RNAi)-treated animals. Several soaking protocols for dsRNA uptake were investigated using the fluorescent dye FITC. Desiccation-enhanced soaking showed the strongest uptake of FITC and resulted in a significant and consistent decrease of mRNA levels of two of the four tested genes (rps-2 and tps-2a), suggesting effective uptake of dsRNA-containing solution by the nematode. These findings suggest that RNAi by desiccation-enhanced soaking is viable in Panagrolaimus sp. DAW1 and provide the first functional genomic approach to investigate freezing tolerance in this non-model organism. RNAi, in conjunction with qPCR, can be used to screen for candidate genes involved in intracellular freezing tolerance in Panagrolaimus sp. DAW1

Molecular analysis of the cold tolerant Antarctic nematode, Panagrolaimus davidi

Isolated and established in culture from the Antarctic in 1988, the nematode Panagrolaimus davidi has proven to be an ideal model for the study of adaptation to the cold. Not only is it the best-documented example of an organism surviving intracellular freezing but it is also able to undergo cryoprotective dehydration. As part of an ongoing effort to develop a molecular understanding of this remarkable organism, we have assembled both a transcriptome and a set of genomic scaffolds. We provide an overview of the transcriptome and a survey of genes involved in temperature stress. We also explore, in silico, the possibility that P. davidi will be susceptible to an environmental RNAi response, important for further functional studies

Thermal reaction norms and the scale of temperature variation: latitudinal vulnerability of intertidal Nacellid limpets to climate change

The thermal reaction norms of 4 closely related intertidal Nacellid limpets, Antarctic (Nacella concinna), New Zealand (Cellana ornata), Australia (C. tramoserica) and Singapore (C. radiata), were compared across environments with different temperature magnitude, variability and predictability, to test their relative vulnerability to different scales of climate warming. Lethal limits were measured alongside a newly developed metric of “duration tenacity”, which was tested at different temperatures to calculate the thermal reaction norm of limpet adductor muscle fatigue. Except in C. tramoserica which had a wide optimum range with two break points, duration tenacity did not follow a typical aerobic capacity curve but was best described by a single break point at an optimum temperature. Thermal reaction norms were shifted to warmer temperatures in warmer environments; the optimum temperature for tenacity (Topt) increased from 1.0°C (N. concinna) to 14.3°C (C. ornata) to 18.0°C (an average for the optimum range of C. tramoserica) to 27.6°C (C. radiata). The temperature limits for duration tenacity of the 4 species were most consistently correlated with both maximum sea surface temperature and summer maximum in situ habitat logger temperature. Tropical C. radiata, which lives in the least variable and most predictable environment, generally had the lowest warming tolerance and thermal safety margin (WT and TSM; respectively the thermal buffer of CTmax and Topt over habitat temperature). However, the two temperate species, C. ornata and C. tramoserica, which live in a variable and seasonally unpredictable microhabitat, had the lowest TSM relative to in situ logger temperature. N. concinna which lives in the most variable, but seasonally predictable microhabitat, generally had the highest TSMs. Intertidal animals live at the highly variable interface between terrestrial and marine biomes and even small changes in the magnitude and predictability of their environment could markedly influence their future distributions

Molecular snapshot of an intracellular freezing event in an Antarctic nematode

The Antarctic nematode, Panagrolaimus sp. DAW1 (formerly called Panagrolaimus davidi), is the best documented example of an organism able to survive intracellular ice formation in all of its compartments. Not only is it able to survive such extreme physiological disruption, but it is able to produce progeny once thawed from such a state. In addition, under slower rates, or less extreme degrees, of cooling, its body remains unfrozen and the vapour pressure difference between the supercooled body fluids and the surrounding ice leads to a process termed cryoprotective dehydration. In contrast to a fairly large body of work in building up our molecular understanding of cryoprotective dehydration, no comparable work has been undertaken on intracellular freezing. This paper describes an experiment subjecting cultures of Panagrolaimus sp. DAW1 to a range of temperatures including a rapid descent to −10 °C, in a medium just prior to, and after, freezing. Through deep sequencing of RNA libraries we have gained a snapshot of which genes are highly abundant when P. sp. DAW1 is undergoing an intracellular freezing event. The onset of freezing correlated with a high production of genes involved in cuticle formation and subsequently, after 24 h in a frozen state, protease production. In addition to the mapping of RNA sequencing, we have focused on a select set of genes arising both from the expression profiles, as well as implicated from other cold tolerance studies, to undertake qPCR. Among the most abundantly represented transcripts in the RNA mapping is the zinc-metalloenzyme, neprilysin, which also shows a particularly strong upregulated signal through qPCR once the nematodes have frozen

Extremophiles in an Antarctic Marine Ecosystem

Recent attempts to explore marine microbial diversity and the global marine microbiome have indicated a large proportion of previously unknown diversity. However, sequencing alone does not tell the whole story, as it relies heavily upon information that is already contained within sequence databases. In addition, microorganisms have been shown to present small-to-large scale biogeographical patterns worldwide, potentially making regional combinations of selection pressures unique. Here, we focus on the extremophile community in the boundary region located between the Polar Front and the Southern Antarctic Circumpolar Current in the Southern Ocean, to explore the potential of metagenomic approaches as a tool for bioprospecting in the search for novel functional activity based on targeted sampling efforts. We assessed the microbial composition and diversity from a region north of the current limit for winter sea ice, north of the Southern Antarctic Circumpolar Front (SACCF) but south of the Polar Front. Although, most of the more frequently encountered sequences were derived from common marine microorganisms, within these dominant groups, we found a proportion of genes related to secondary metabolism of potential interest in bioprospecting. Extremophiles were rare by comparison but belonged to a range of genera. Hence, they represented interesting targets from which to identify rare or novel functions. Ultimately, future shifts in environmental conditions favoring more cosmopolitan groups could have an unpredictable effect on microbial diversity and function in the Southern Ocean, perhaps excluding the rarer extremophiles

Characterisation of the mantle transcriptome and biomineralisation genes in the blunt-gaper clam, Mya truncata

Members of the Myidae family are ecologically and economically important, but there is currently very little molecular data on these species. The present study sequenced and assembled the mantle transcriptome of Mya truncata from the North West coast of Scotland and identified candidate biomineralisation genes. RNA-Seq reads were assembled to create 20,106 contigs in a de novo transciptome, 18.81% of which were assigned putative functions using BLAST sequence similarity searching (cuttoff E-value 1E − 10). The most highly expressed genes were compared to the Antarctic clam (Laternula elliptica) and showed that many of the dominant biological functions (muscle contraction, energy production, biomineralisation) in the mantle were conserved. There were however, differences in the constitutive expression of heat shock proteins, which were possibly due to the M. truncata sampling location being at a relatively low latitude, and hence relatively warm, in terms of the global distribution of the species. Phylogenetic analyses of the Tyrosinase proteins from M. truncata showed a gene expansion which was absent in L. elliptica. The tissue distribution expression patterns of putative biomineralisation genes were investigated using quantitative PCR, all genes showed a mantle specific expression pattern supporting their hypothesised role in shell secretion. The present study provides some preliminary insights into how clams from different environments – temperate versus polar – build their shells. In addition, the transcriptome data provides a valuable resource for future comparative studies investigating biomineralisation

Molecular mechanisms underpinning transgenerational plasticity in the green sea urchin \u3cem\u3ePsammechinus miliaris\u3c/em\u3e

The pre-conditioning of adult marine invertebrates to altered conditions, such as low pH, can significantly impact offspring outcomes, a process which is often referred to as transgenerational plasticity (TGP). This study describes for the first time, the gene expression profiles associated with TGP in the green sea urchin Psammechinus miliaris and evaluates the transcriptional contribution to larval resilience. RNA-Seq was used to determine how the expression profiles of larvae spawned into low pH from pre-acclimated adults differed to those of larvae produced from adults cultured under ambient pH. The main findings demonstrated that adult conditioning to low pH critically pre-loads the embryonic transcriptional pool with antioxidants to prepare the larvae for the “new” conditions. In addition, the classic cellular stress response, measured via the production of heat shock proteins (the heat shock response (HSR)), was separately evaluated. None of the early stage larvae either spawned in low pH (produced from both ambient and pre-acclimated adults) or subjected to a separate heat shock experiment were able to activate the full HSR as measured in adults, but the capacity to mount an HSR increased as development proceeded. This compromised ability clearly contributes to the vulnerability of early stage larvae to acute environmental challenge